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Abstract 

Neutrino emission caused by Cooper pairing of baryons in neutron 
stars is recalculated by accurately taking into account for conservation of 
the vector weak current. The vector current contribution to the neutrino 
emissivity is found to be several orders of magnitude smaller than that 
obtained before by different authors. Therefore, the neutrino energy losses 
due to singlet-state pairing of baryons can in practice be neglected in 
simulations of neutron star cooling. This makes negligible the neutrino 
radiation from pairing of protons or hyperons. The neutrino radiation 
from triplet pairing takes place through axial weak currents. For these 
states, when the total momentum projection is rtij = 0, the vanishing 
of the vector weak current contribution results in the suppression of the 
neutrino energy losses by about 25%. The neutrino emissivity due to 
triplet pairing with \rrij \ = 2 is suppressed by about a factor of 3, caused by 
the collective contribution of spin-density fluctuations in the condensate. 

1 Introduction 

Thermal excitations in superfluid baryon matter of neutron stars, in the form 
of broken Cooper pairs, can recombine into the condensate by emitting neu- 
trino pairs via neutral weak currents. The neutrino energy losses caused by the 
pairing of neutrons in a singlet state was first calculated by Flowers et al. [1 
and reproduced by other authors 0.0 Yakovlev et al. p] have also consid- 
ered the neutrino energy losses in the case of neutron pairing in a triplet state. 
Jaikumar & Prakash ^] have generalized this mechanism to superconducting 
quark matter. The neutrino energy losses due to pairing of hyperons pj, are 
also discussed in the literature as possible cooling mechanisms for superdense 
baryonic matter in neutron stars. Nowadays, these ideas are widely accepted 
and used in numerical simulations of neutron star evolution [Jj, jH], 0. Never- 
theless, the existing theory of neutrino radiation from Cooper pairing of baryons 
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(quarks) in neutron stars rises some questions because, as we show below, it is 
inconsistent with the hypothesis of conservation of the vector current in weak 
interactions. 

Let us recall shortly the main steps in the above calculations. The low- 
energy Lagrangian of the weak interaction may be described by a point-like 
current-current approach. For interactions mediated by neutral weak currents, 
it can be written as 1 

Gf 

Here Gf is the Fermi coupling constant, and the neutrino weak current is given 
by — (1 — 75) v. The vacuum weak current of the fermion (baryon or 
quark) is of the form J M = ^ (CyJ^ — CaJ^Js) tp, where, V represents the 
fermion field, and the weak vertex includes the vector and axial-vector terms 
with the corresponding coupling constants Cy and G'a- 

Since relativistic calculations are more complicated and less transparent, 
we consider the nonrelativistic case, typical for superfluid baryon matter in 
neutron stars. Then, the nonrelativistic limits for the baryon operators are 
■ipB^ipB — ► ^b^b, i>Blil^i>B — > &g&i&B, all others being zero. Here ^b is 
the second-quantized nonrelativistic spinor wave function, and &i are the Pauli 
matrices. 

The process is kinematically allowed due to the existence of a superfluid 
energy gap A, which admits the transition with time-like momentum transfer 
K = (iu, k), as required by the final neutrino pair. We have to > 2 A and w > k . 
The emissivity for neutrino pairs due to recombination of baryon quasi-particle 
excitations in the singlet state is then 2 

G F \ 2 2 f d 3 pd 3 p' fd 3 k 1 d 3 k 2 2 2 

x6(p + p' -k) 6 (e p + e p , -lu). (2) 

Here u> = u>i + u>2 and k = ki + k2 are the energy and momentum carried out by 
the freely escaping neutrino pair, and / (e p ) is the Fermi distribution function 
of the quasi-particles with energy e p , as given by Eq. ifTTjl . The neutrino matrix 
element is of the standard form \M. V \ = 8 (uilu 2 — kik 2 ). Therefore, the rele- 
vant input for this calculation is the recombination matrix element between the 
baryon state, which has a pair of quasi-particle excitations of momentum-spin 
labels (p, up; p', down), and the same state but with these excitations restored 
to the condensate. To the leading (zero) order in k <C pf, this matrix element 
is usually estimated as 

\M B \ 2 - ^ (3) 

v 

yielding the following neutrino energy losses at temperature T < T c . 

n„\ AG%p F M*Cl Kr m7 2 f°° z*dx 
Qfrs ( So) = — 3-5 -N„T'y / 

1 In what follows we use the Standard Model of weak interactions, the system of units h = 
c = 1 and the Boltzmann constant ks = 1. The fine-structure constant is a = e 2 /4ir = 1/137. 

2 In the case of 1 So pairing, the total spin of the Cooper pair is zero, and the axial contri- 
bution vanishes in the nonrelativistic limit. 
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where M* is the effective nucleon mass, y — A/T, z — 1 x 1 + y 2 , and M v = 3 is 
the number of neutrino flavors. T c is the critical temperature for baryon pairing. 

The naive estimate 10 is inconsistent with the hypothesis of conservation of 
the vector current in weak interactions. Indeed, a longitudinal vector current 
of quasi-particles consisting only on a temporal component can not satisfy the 
continuity equation. 

The question of conservation of the longitudinal vector current caused by 
recombination of quasi-particles has been discussed by many people ^01 m con- 
nection with the gauge invariance of the Bardeen-Cooper-Schrieffer theory of 
superconductivity. It was realized that the current conservation could be re- 
stored if the interaction among quasi-particles is incorporated in the coupling 
vertex to the same degree of approximation as the self-energy effect is included 
in the quasi-particle. It has been also pointed out that there is significant differ- 
ence between the transverse and longitudinal current operators in their matrix 
elements. Namely, there exist collective excited states of quasi-particle pairs 
[TT] . [12] which can be excited only by the longitudinal current. As a result, 
the spatial part of the longitudinal current does not vanish, and cancels the 
temporal part. 

In the present paper we recalculate the neutrino energy losses with allowance 
for conservation of the weak vector current. The paper is organized as follows. 
In section 2 we shortly discuss the Nambu-Gorkov formalism as a convenient 
description of the particle-hole excitations in the system with pairing. The 
wave functions of the quasi-particle excitations are obtained in section 3. In 
section 4 we derive an effective vertex conserving the vector weak current for 
neutral and charged baryons. In section 5 we use the quasi-particle states for 
the calculation of the matrix element of the effective vector weak current, and 
the neutrino energy losses in the vector channel. The neutrino energy losses 
in the axial channel are calculated in section 6. A discussion of the obtained 
results and our main conclusions are presented in section 7. 



2 Formalism 

In the Nambu-Gorkov formalism, the quasi-particle fields are represented by 
two-component objects 



i>l (-p) 

if>i (p) is the the quasi-particle component of the excitation with momentum p 
and spin a, and ?/4 (~P) ls the hole component of the same excitation, which can 
be interpreted as the absence of a particle with momentum — p and spin — er. 
The two-component fields (JjJ obey the standard fermion commutation relations 

With the aid of the 2x2 Pauli matrices 

fl = ( 1 ) ' fa = ( i o* ) ' h = ( -1 ) (6) 

operating in the particle-hole space, the Hamiltonian of the system of quasi- 
particles can be recast as |12| 

H = Hq + If i , 
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where 



#o = £>J A (i>))* P (7) 



is the BCS reduced Hamiltonian, where the nonrelativistic energy is measured 
relatively to the Fermi level 

t - P - 

Sp - 2M * Mj 

M* is the effective mass of the quasi-particle and /i is the Fermi energy. The 
quasi-particle self-energy, of the form 

is a 2 x 2 matrix in the Nambu-Gorkov space, and a 2 x 2 matrix in the spin space, 
which depends on the orientation of the quasi-particle momentum p = p/ |p| 
(see e.g. [T5]). 

The residual interaction among quasi-particles is given by the following 
Hamiltonian 

#1 = \ E V PP' (?) (*P+ 9 f 3*p) (^-/3*p/ 

p,p' 

As follows from the Hamiltonian JJJ ; the inverse of the quasi-particle prop- 
agator has the simple form: 

^ 1 b)= P o-e P T 3 -A(p) (9) 



3 Quasi-particle states 

Near the Fermi-surface, the imaginary part of the quasi-particle energy is small. 
Therefore, to the extent that the single-particle picture makes some physical 
sense, we will describe the quasi-particles with the aid of wave-functions. The 
states of quasi-particles obey the equation 

Q~ 1 V P = 0. (10) 

Writing ^ p in the form 

with a and being two-component spinors, we arrive to the following set of 
matrix equations 

(po-Si P )a-Z)(3 = 
&a-(p + Q(3 = 

From the second equation we have 

a = ( P o + Hp) (D^y 1 13 (12) 
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By substituting this into the first equation and multiplying from the left, we 
obtain 

(po - £p - &D) = (13) 
This equation has nontrivial solutions only if 



det 



p 2 - £ - = 0. (14) 



In the interesting cases of singlet- and triplet-state pairing we are going to 
consider, the gap matrix D is proportional to the unitary matrix T (p) 



£(p) = A p T(p) (15) 
a scalar function Q^l, so that 



with Ap 



f tf = 1 (16) 
Eq. ((Hj) then gives us the eigenvalues 



po = ±6 P , e p = ^/Q + A? (17) 

with 

A? = iTrftD (18) 

Using the fact that Eq. ((13)) is diagonal, we may choose two independent 
solutions for /? as the ordinary spinors: 

= ( I ), /3 2 ± = x;=(5) (19) 

Here, the upper sign corresponds to the positive- or negative-frequency solution, 
in accordance with Eq. ifTTjl . The a-component of these solutions is to be found 
from Eq. ltl2)l . 

From Eq. 11511 . we obtain 

(z^^A^f (p) 

Substitution this into Eq. Ill2t gives 

af = A7 1 (±e p +^)f/3 i (20) 

Making use of Eqs. l(T9|) . ll2f))) we can find the normalized wave functions. 
The positive-frequency states can be written as 



q, = J_( ^(tp+^TXa \ e ipr-i6 p t 



pa /N 



(21) 



The factor in Eq. Illl|l is to be found from the normalization condition 

^p^p = 1. A direct evaluation making use of Eq. ljTfi)l yields 



1 = / e P - £p 
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Incorporating this factor into Eq. (12H we obtain the positive-frequency eigenstates 



* p >a = ( U P TXa ) e lpr - lt ^ (22) 



VpXa 



with 





(23) 



The negative-frequency wave-functions can be obtained in the same way. The 
state with momentum p and spin-label a has the form 

¥- p ,-a = ( ~ v P Tx - a ) e - ipr+ie ^ (24) 

This solution is connected to the hole state by the particle-antiparticle conju- 
gation 

C : * c = C** = f 2 *t. 

which changes quasi-particles of energy-momentum (jjq, p) into holes of energy- 
momentum (— pq, — p), or interchanges up-spin and down-spin particles. 

It is necessary to note that, in general, ^p,^ is not an eigenstate of the 
spin projection, because the upper component oc T\a is a linear combination of 
different spin states. The label a indicates only the spin function of the lower 
component, and serves for identification of the quasi-particle states. 



4 Effective vertex for quasi-particles 

4.1 Neutral baryons 

The components of the bare vertex 

7* = ( f 3 V = ' inn ■ (25) 

I j^P if ju = i = 1,2,3 

are 2x2 matrices in the Nambu-Gorkov space. As already mentioned, the longi- 
tudinal current corresponding to the bare vertex does not satisfy the continuity 
equation. To restore the current conservation, one must consider the modifica- 
tion of the vertex 7^ to the same order as the modification of the propagator 
is done. The relation between the modified vertex T^ and the quasi-particle 
propagator (O is given by the Ward identity JBj 

KpT* ( P ',p) = faG- 1 (p) - G- 1 ip') f 3 , (26) 

where K = (u, k) is the transferred momentum. The plane wave solutions 

!i Pi(1 exp (ipr — if.pt) , u* p t a , exp (— ipr + ie p t) 

for W and \P + obey the equations G -1 (p)u PlCt — 0, and up/^Q^ 1 (p') = 0. 
Therefore the Ward identity implies conservation of the vector current on the 
energy shell of the quasi-particles. 
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Figure 1: Dyson equation, in the ladder approximation, for the effective vertex. 
Solid lines stand for dressed particles, and the dashed line represents the Z-boson 
field. The wavy line accounts for the residual interaction among quasi-particles. 



Following the prescriptions of quantum electrodynamics, an approximation 
which satisfies the Ward identity (and hence the continuity equation) is the sum 
of ladder diagrams satisfying the Dyson equation shown in Fig. 1. 

In these diagrams, all solid lines are "dressed" particles, and the dashed line 
is the Z-boson field. The "dressed" particles interact with the same primary 
interaction V pp > which produces the self-energy of the quasi-particle. Thus, for 
K = p — p' , the corrected vertex satisfies the integral equation 

T»(p-K,p) = r(p~K,p) (27) 

+* J r 3 Q (p' - K) I* {p' - K,p') Q [p') t 3 V pp , 

In the limit K = (uj, 0), the Ward identity gives 

r°(p-iC,p)=f 3 --f 3 A(p) 

with 



ui 



The poles of the vertex function correspond to collective eigen-modes of the 
system. Therefore, the pole which appears at w — > 0, k — implies the existence 
of a collective mode, which plays an important role in the conservation of the 
vector current. The corresponding nonperturbative solution to Eq. 112 7JI has 
been found by Nambu [12] (see also ^lj). In our notation, it reads 3 

r (p - K, p) = f 3 - 2f 3 A (p) (29) 

The poles in this vertex correspond to the collective motion of the condensate, 
with the dispersion relation ui — ak, where a 2 = Vp/3. 

The effective vertex satisfies the Ward identity 12611 , and thus the continuity 
equation on the energy shell 

uTn - kT ~ (31) 

3 To obtain the weak vector current, this vertex should be multiplied by the weak coupling 
constant Cy- 
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Figure 2: Dyson equation for the vertex correction for charged quasi-particles 
(compare with Fig. 1). The shaded areas represent the modified effective vertex, 
and the wavy line stands for the Coulomb interaction. 

Indeed, by making use of ui = e p + e p > and k = p + p' with k <C p ~ pp, 
the left-hand side of this equation can be reduced to the right-hand side of the 
Ward identity Iplj) . 

4.2 Charged baryons 

Consider now the case when quasi-particles carry an electric charge. Including 
the long-range Coulomb interaction Vc (k) — e 2 /k 2 implies that the vertex part 
is multiplied by a string of closed loops, which represents the polarization of 
the surrounding medium. In this case, the new vertex can be found as 
the solution of the Dyson equation, according to the diagram of Fig. 2. or, 
analytically 



(p-K,p) 



(p -K,p)- r - K,p) V c (fc) 
Xl I 7^ Tr h<? (P' - K ) ^ (P' - K >P') (P')} ( 32 ) 




The integral 



ft " (k) 




1 Tr t 3 Q (p' - K) f" {p' - K,p') Q (p') 



(33) 



in the r.h.s. of Eq. 1(32)1 represents the n° M component of the polarization tensor 
shown by the loop in Fig. 2. Introducing also the notation 




Eq. l|.'S2ll can be recast as follows 



rr> (K) = n ^ 1 (k) - v c (k) n 00 (K) fi ^ (k) . 



We then obtain 



ft ' 1 {K) 



n ^ 1 (k) 



i + v c (k) n 00 {K) ■ 



Combining this with Eqs. (l3~2"l) . (HT^ft yields 



f" (p - K,p) = (p - K,p) 1 



V c (k) n ^ (K) 



) 



1 + V c (k) n 00 (K) 
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In particular, for T° we arrive to 4 



,P) l + V c {k)Il™(K)- 



(34) 



The polarization function II 00 (K) can be readily calculated with the help 
of r° given by Eq. (1201 . By neglecting the small dependence of the energy gap 
on the transferred momentum k, we have 



n °° ( K ) = * I 44 Tr [f 3 Q (p' - K) f 3 Q (p')} 

d 4 p' 



2lo 



Tr 



u 2 -a 2 k 2 "J (2tt) 4 
Here the quasi-particle propagator follows from Eq. Q : 



r 3 Q(p' -K)r 3 A(p')Q(p') 



G (P) = 2 l _ e 2 (P0 + CpT3 + A (p)) 



(35) 



Performing the trace over particle-hole and spin indices, and integration over 
Po, results in the following expression 



n uu (u, k) 



(2tt) 3 e P (w + e p + £p-k) (w - e p -k + e P ) 

We are interested in the regime defined by 

k < u), io > 2e p » e p — e p _k — fcVp 
In this case we obtain, after some simplifications 5 



2A? 



2d 3 p uj€ p ^ k - A P - £ P gp-k - ep_ k - a A k 2 [u 2 - a 2 k 

(2tt) 3 e P -k (w - e p - e p _k) (w + e p - e p _ k ) 

2rf 3 p - W £ P - A| - e p - C P £p-k - a 2 fc 2 (c 2 - a 2 ^)" 1 2A 2 



n uu (-,k) = — 



1 

8^ 3 



2 1,2 



a k 



a 2 k 2 



d 3 p 



-a 2 fc 2 /4) 



fc 2 / d 3 p 



ip (uj 2 



with a 2 = V|/3. 

Inserting this into Eq. 1)3411 gives 

f°(p-K,p) 



3M 2 e 2 



a 2 fc 2 /4 



T°(p-K,p) 



(36) 



4 The solution to the equation I + Vq (fc) II 00 (A") = determines the new dispersion law 
oj = u> (fc) for the collective excitations, which represents plasma waves |12) . 

6 An analogous expression has been obtained in 112J for the case of singlet pairing of elec- 
trons. 
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with 



a 2 k 2 



d 3 p 



+ / d 3 p 



.,,2 



2 fc 2 /4) 



a 2 fc 2 /4 



(37) 



5 Energy losses in the vector channel 
5.1 Neutral baryons 

Having at hand the effective vertex and the wave functions of quasi-particles 
and holes, we can evaluate the matrix element of the vector weak current. In 
the particle-hole picture, the creation and recombination of two quasi-particles 
is described by the off-diagonal matrix elements of the Hamiltonian, which cor- 
responds to quasi-particle transitions into a hole (and a correlated pair). Thus, 
we calculate the matrix element of the current between the initial (positive- 
frequency) state of a quasi-particle with momentum p and the final (negative- 
frequency) state with the same momentum p. 

M li = {*l. Pt _„\T lt \* Vtir ) 

Let us consider separately the contributions from the bare vertex, given by 
the first terms in Eq. p9|l. and the collective part, given by the second term, 
so that M» = M^ vc + M c ° l \ 

Making use of the wave functions described by Eqs. 12211 . Il24t and using the 
identity 11 fill , for /x — 0, we find 

Mo arc = - K V + W ) ^ - — , k<^p~p F (38) 

with e p + e p i — uj and p + p' = k. 

The velocity of the collective mode a 2 = Vp /3 is small in the nonrelativistic 
system. Therefore, we expand the collective contribution in this parameter to 
obtain 

The contribution of the bare vertex .Mo are reproduces the matrix element 
<EJ derived by Flowers et al. [I] and Yakovlev et al. [jjj. However, the collective 
correction modifies this crucially. In the sum of the two contributions, the 
leading terms mutually cancel, yielding the matrix element 

Mo = + M c oU ~ ^V?4 — 

which is at least ~ Vp times smaller than the bare result. 

The above cancellation of the contribution of "bare" vertex is not acciden- 
tal. This part of the matrix element, which remains finite at zero transferred 
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momentum, cannot be restored with taking into account of more complicated 
corrections to the vertex part. It is well-known that, for any process involving 
interactions, the perturbation diagrams can be grouped into gauge invariant 
subsets, such that the invariance is maintained by each subset taken as a whole. 
If any new subset of the diagrams is added to the vertex part, the same should 
also be incorporated into the quasi-particle self-energy. If this is carefully done, 
for k — > 0, one must obtain the exact cancellation of the contribution of the 
"bare" vertex. Indeed, ^S> + (r, t) (r, t) represents the operator of baryon charge 
density. If the ground state |0) is an eigenstate of the baryon charge operator 
for the system with N particles, and \n > is a state for the system with N + 2 
particles, the matrix element 

J (0| *+ (r, t) * s (r, t) \n) e lkr - luJt d 3 rdt 

should vanish for k — > 0, because the above expression is an off-diagonal matrix 
element of the total charge operator Q. This implies that, for k <C pf , the 
matrix element should be proportional to some power of k. 

The spatial component of the longitudinal (with respect to k) component 
of the matrix element can be obtained from Eq. 13 111 . Since kT = (u>/k) Tq we 
have 

■M-W — o V F • 

3 u> e p 

In the above, k = k/k is a unit vector directed along the transferred momentum. 

Since the collective interaction modifies only the longitudinal part of the 
vertex, the transverse part of the matrix element can be evaluated directly from 
the bare vertex <25li . This yields 

M±_ ~ (VpU p , - u p v p ,) ~ _Iy|^£ (kp) p ± . 

The rate of the process is proportional to the square of the matrix element. 
This means that the vector current contribution to the neutrino energy losses is 
Vp times smaller than estimated before. The corresponding neutrino emissivity 
in the vector channel can be evaluated with the aid of Fermi's golden rule: 

2^T^f wTr M * M " S (p + p' - k ) 6 ( £ p + £ p' - w ) • 

One can simplify this equation by inserting J d 4 K 5^ (K — k\ — fca) = 1- Then, 
the phase-space integrals for neutrinos are readily done with the aid of Lenard's 
formula 

1*^(4) {K _ kl _ k2)Tr(m 

ZUJi ZLO2 

Aqr 

= — (Kp,K„ - K 2 g, v ) 6 (K 2 ) e H , 
where Q(x) is the Heaviside step function. 
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For k <C Pf we obtain 

QV = ~3~{^k) ^nf K J d3pf2{£p) Jo ^ Jo 
x ((^M") 2 - if iM^) <J (2e p - u) . 

The next integrations are trivial. We get 



-V*G 2 F Cl PF M*T 7 y 2 j 



This is to be compared with Eq. 0|, We see that 

Qv 148 



Qfrs (iS°) 2835 



V, 
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i.e. the neutrino radiation via the vector weak currents in the nonrelativistic 
system (Vf <C 1) is suppressed by several orders of magnitude with respect to 
the predictions of Flowers et al. pp. 



5.2 Charged baryons 

If the paired baryons carry an electric charge, the effective vector vertex is given 
by Eqs. ffi6jl . I|29|) and l|37jl. We are interested in the case u> > 2Ap and k < u>. 
Since o« 1 and p ~ pf, the second integral in Eq. I|37|) may be dropped. By 
neglecting also the small contributions from a 2 k 2 <C e 2 ,w 2 we get 

X( ] W e^ p ( 2 ^) 3 ^M*J Itt J 
where n is the number of baryons per unit volume. Since 

J 4tt J e 3 J 4tt J a e 2 ^/ £ 2 _ A 2 

J 4^ ^ x 2 V^T 

we obtain 

This agrees with the plasma frequency for a free gas of charged particles. 

The energy exchange in the medium goes naturally as the temperature scale. 
Therefore, the energy transferred to the radiated neutrino-pair is u> ~ T < 
T c , while the plasma frequency oj p is typically much larger than the critical 
temperature for Cooper pairing. For instance, for a number density n of the 
order of the nuclear saturation density no — 0.17 fm and the effective mass of 
the baryon M * of the order of the bare nucleon mass, we obtain uj p ~ 10 MeV, 
while the critical temperature for baryon pairing is about 1 MeV or less. Under 
these conditions, we obtain 
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f° (p - k, p ) ~ — r° (p - jf, P ) ~ -^r° (p - *r,p) 

Thus, in superconductors, the vector current contribution to the neutrino radia- 
tion is suppressed additionally by a factor (T^/oj^) : this is the plasma screening 
effect. The total suppression factor, due to both the current conservation and 
the plasma effects, is of the order 

(T?H) V F 4 < 10" 6 . 



6 Energy losses in the axial channel 

The neutrino energy losses in the axial channel can also be obtained with the aid 
of Fermi's golden rule. After integration over the phase space of participating 
particles, this yields (see details in 



15tt 5 8tt J J ( e * + \y 



Here M v — 3 is the number of neutrino flavors, z — J x 2 + A?/T 2 , and the 
outer integration is performed over the orientations of the nucleon momentum 
p. The matrix elements of the axial current 



come into this expression under the combination 

I = Ixx+ Iyy + Izz (41) 

with 

hi = Y^ A i it>,ri)A*{p,r,') 

nv' 

The collective mode considered in previous sections represents density os- 
cillations of the condensate, and does not influence the axial weak vertex of a 
quasi-particle, if the average spin projection of the bound pair is zero. This 
happens in the singlet-state pairing and in triplet pairing with total momentum 
projection rrij = 0. In these cases, the axial vertex may be taken in its bare 
form A\ are = Sj, where the spin operator 

4), (42) 

with 6i = (<t x , a y , a z ) being the Pauli spin matrices, is defined according to the 
relation 

where 



VpX-v 
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is the eigenstate with spin projection s z = rj/2 = ±1/2 along the Z-axis. 

A direct evaluation of the matrix elements of the bare axial vertex 1421 for 
a quasi-particle transition into a hole by making use of the wave functions l!22t , 
i|24l gives I — in the case of singlet pairing. For triplet pairing with mj = 
we obtain 

A? A? 

Ixx = Iyy = (1 - COS 2(f) , I zz = 2—f 

which, after averaging over the azimuthal angle ip of the quasi-particle momen- 
tum p, gives 

A? 

J = 4^. (43) 
e P 

This confirms the corresponding result obtained in pj]. 

In the case of a triplet-state pairing with maximal momentum projection 
\nij\ = 2, we obtain a different result. In this case, the spin projection of the 
Cooper pair is conserved, and collective density oscillations of the condensate 
are accompanied by oscillations of the spin density, i.e. by the axial current. 
The corresponding gap matrix of the quasi-particles has the form D — A p T (p) 
with [H] 

/ — e l f \ 
T(p)=l e _ lip J = —cr z cos <p — i sin ip. 

Therefore, the operator S z commutes with the quasi-particle Hamiltonian Hq = 
£,pT3 + A(p), where A is given by Eq. JSJ. This modifies the z-component of 
the axial vertex, which should be found from the equation 



A z (p - K,p) = r 3 a z + i f r 3 Q {p' - K) A z (p' - K,p r ) Q (p') r 3 V pp , 

J [2tt) 



After multiplication of this equation by a z from the right-hand side, and intro- 
ducing the new function 

B = A z (p-K,p)a z (44) 
we arrive to the following equation 

B{p-K,p)=T 3 + i J 0^ f 3 Q (p' - K) B (p' - K,p') a z Q (p') a z f 3 V pp ,. 

By using the fact that a z commutes with Q{p'), one can reduce this equation 
to the form 



B(p- K,p) = r 3 +i f t 3 Q (P — K)B {p' - K,p') Q {p') r 3 V vv , 

J (27T) 



which is identical to Eq. i|2Tt for p — 0. Then, the solution has a form analogous 
to Eq. (EH!)- We obtain 

B(p-K,p)=T 3 -2T 3 k(p) 



uj 2 — a 2 k 2 
From Eq. lUijl we find 

A z (p - K,p)= r 3 a z - 2f 3 A (p) a 2 



lo 2 — a 2 k 2 ' 
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or, in an explicit form: 

o z \ _ ( & z f \ LoA p 



A z {p-K,p) 



-a, J \ -«r z Tt J uj 2 -a 2 k 2 ' 
The transverse components of the axial vertex are not modified 

^ = ^=(^0 -* x )> A y = T ^y={o -a v 

Direct evaluation of the matrix elements for k <C pf gives, in this case: 
A x (p ! r?) = -^(l-e 2 w)^, 

^ ( P , „) = (i _ ^) («y ti , - 6 Un ) . 

Due to the contribution of the collective mode, the matrix element A z is Vp 
times smaller than A x and A y and, thus, can be neglected. After averaging over 
the azimuthal angle cp of the quasi-particle momentum, we obtain 

I~I XX + Iyy = 2-£ (45) 

Thus, in the case of 3 Pi pairing with total momentum projection \m,j\ = 2, 
the value of /, given by Eq. ljl5]). is twice smaller, and the neutrino energy 
losses are proportionally suppressed. 



7 Discussion and conclusions 

We have considered the problem of conservation of the vector weak current 
in the theory of neutrino-pair radiation from Cooper pairing in neutron stars. 
The correction to the vector weak vertex is calculated within the same order 
of approximation as the quasi-particle propagator is modified by the pairing 
interaction in the system. This correction restores the conservation of the vector 
weak current in the quasi-particle transition into the paired state. As a result, 
in the nonrelativistic baryon system, the matrix element of the vector current is 
Vp times smaller than previous estimations. This means that the vector weak 
current contribution to neutrino radiation caused by Cooper paring is Vp times 
smaller than it was thought before. The vector weak current contribution from 
pairing of charged baryons is suppressed additionally by a factor ~ (T c 2 jujfy 
due to plasma screening. The total suppression factor due to both the current 
conservation and the plasma effects is of the order 

(T 2 J^l) V F 4 < 10" 6 . 

In papers ^2], |17j . a special case of proton pairing has been considered 
which accounts for electron-proton correlations, resulting in an increasing of the 



15 



neutrino emissivity. This, however, takes place only because of a very small weak 
vector coupling constant of the proton with respect to that for the electron. In 
our calculations, incorporation of the electron-proton correlations would result 
in the replacement of the vector weak coupling constant of the proton by a 
form-factor, which is proportional (to the leading order) to the electron coupling 
constant. This could enlarge the effective proton vector weak coupling to almost 
the same order of magnitude as that for neutrons. However, the cancellation of 
the temporal component of the transition current due to current conservation 
is much more stronger, and makes negligible the neutrino emission from the 
proton pairing even with taking into account the electron-proton correlations. 
Thus, the neutrino energy losses due to singlet-state pairing of baryons can, 
in practice, be neglected in simulations of neutron star cooling. This makes 
unimportant the neutrino radiation from pairing of protons or hyperons. 

The neutrino radiation from triplet pairing occurs through the axial weak 
currents. In the case of 3 P2 pairing with total momentum projection rrij = 0, 
the corresponding neutrino emissivity is given by 

Q K = 0) = c f G2 f/ 5 M * WF u (46) 

107T 

with 



1 f in 2 [°° , {x z +y 



2 \ 



p t = TZ / dSl y 2 / dx \ °- ' s2 . (47) 



This formula can be obtained directly from the corresponding expression sug- 
gested in p] by omitting the contribution proportional to the vector weak cou- 
pling constant C v . As a result, the corresponding neutrino energy losses are 
suppressed by about 25% 

Q (rrij = 0) _ 2C\ _ 3.17 



<2ykl {rrij = 0) C v + 2C\ 4.17 

with respect to those calculated in 0. 

In the case of 3 P2 pairing with total momentum projection \rrij\ — 2, the 
neutrino energy losses in the axial channel are additionally suppressed due to 
the collective contribution of the spin-density fluctuations in the condensate: 

Q (Kl = 2) = C\ ^ G \ VF f M V T" F t . (48) 

Together with the vanishing of the vector contribution, the total suppression of 
the neutrino emissivity with respect to that obtained in [Hj is about 3 times 

Q(|m 3 | = 2) _ C\ __ 1.58 

Qykl(K I = 2) ~ ci + 2c\ ~ Zlr 

In this paper, we have considered the neutrino radiation from Cooper pair- 
ing in nonrelativistic baryon matter. It is clear from the above consideration, 
however, that conservation of the vector weak current has to be restored also 
in the theory of neutrino radiation from pairing of relativistic quarks |3j. This 
will be done in a future work. 
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